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Nearly half of the exoplanets found within binary star systems reside in very 
wide binaries with average stellar separations beyond 1,000 AU^ (1 AU being 
the Earth-Sun distance), yet the influence of such distant binary companions on 
planetary evolution remains largely unstudied. Unlike their tighter counterparts, 
the stellar orbits of wide binaries continually change under the influence of the 
Galactic tide and impulses from other passing stars. Here we report numerical 
simulations demonstrating that the variable nature of wide binary star orbits 
dramatically reshapes the planetary systems they host, typically Gyrs after for- 
mation. Contrary to previous understanding, wide binary companions may often 
strongly perturb planetary systems, triggering planetary ejections and exciting 
orbital eccentricities of surviving planets. Indeed, observed exoplanet eccen- 
tricities oflfer evidence of this; giant exoplanet orbits within wide binaries are 
statistically more eccentric than those around isolated stars. Both eccentricity 



distributions are well-reproduced when we assume isolated stars and wide bina- 
ries host similar planetary systems whose outermost giant planets are scattered 
beyond ~10 AU from their parent stars via early internal instabilities. Con- 
sequently, our results suggest that although wide binaries eventually truncate 
their planetary systems, most isolated giant exoplanet systems harbor additional 
distant, still undetected planets. 

Unlike binaries with separations below ~10^ AU, very wide binary stars are only weakly 
bound by self-gravity, leaving them susceptible to outside perturbations. As a result, the 
Milky Way's tide and impulses from other passing stars strongly perturb wide binary orbits 
These perturbations, which are fairly independent of the orbiting object's mass, are 
also known to dramatically affect the dynamics of solar system comets with similar orbital 
distances''' ^. Galactic perturbations drive a psuedo-random walk in these comets' pericenters 
(or closest approach distances to the central body)0"^. The same will occur in wide binary 
orbits. Thus, even if a very wide binary's initial pericenter is quite large, it will inevitably 
become very small at some point if it remains gravitationally bound and evolves long enough. 
Such low pericenter phases will produce close stellar passages between binary members, with 
potentially devastating consequences for planetary systems in these binarie^l^. Counterin- 
tuitively, we therefore suspect that wide binary companions could more dramatically affect 
planetary system evolution than tight binaries. 

To investigate this scenario, we use the MERCURY simulation package to perform 2600 
simulations modeling the orbital evolution of our Sun's four giant planets (on their current 
orbits) in the presence of a very wide binary companioiP. These simulations are listed as set 
"A" in Table 1, which briefly summarizes our different simulations' initial conditions (supple- 
mentary information provides details). An example simulation is shown in Figure 1. Initially, 
the binary companion has no effect on the planets' dynamics since its starting pericenter (g) 
is ~3,000 AU. However, after 1 Gyr of evolution, galactic perturbations drive the binary 
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pericenter near 100 AU, exciting the eccentricities of Neptune and Uranus. Once again at 
3.5 Gyrs, the binary passes through another low pericenter phase, this time triggering the 
ejection of Uranus. Finally at 7.2 Gyrs, the binary makes a final excursion to low g, causing 
Neptune's ejection. 

Such behavior is not unusual. Depending on the binary's mass and semimajor axis (mean 
separation, or a,,,). Figure 2a demonstrates that ~30-60% of planetary systems in simulation 
set A experience instabilities causing one or more planetary ejections after 10 Gyrs (the 
approximate age of our galaxy's thin disk). Even though binaries with smaller semimajor 
axes are less affected by galactic perturbations. Figure 2a shows the influence of binary 
semimajor axis on planetary instability rates is weak. This is because tighter binaries make 
pericenter passages at a higher frequency. In addition, when they reach low-g phases they 
remain stuck there for a much longer time than wider binaries. As Figure 2b shows, both of 
these effects cause tighter binaries to become lethal at a much larger pericenter, offsetting the 
Galaxy's diminished influence. The large majority of binary-triggered instabilities are very 
delayed. For binaries with a* > 2, 000 AU, Figure 2c shows that well over 90% of instabilities 
occur after at least 100 Myrs of evolution, well after planet formation is complete. For tighter 
binaries, many more begin in orbits that destabilize the planets nearly instantly. 

While planets are believed to form on nearly circular orbit^, most known giant plan- 
ets (msini > 1 Mjup) have significant non-zero orbital eccentricities (eccentricities of less 
massive planets are known to be colder)'^. This observed eccentricity distribution can be 
reproduced remarkably well when systems of circularly orbiting planets undergo internal 
dynamical instabilities causing planet-planet scattering events that eject some planets and 
excite the survivors' eccentricitie^ For planetary systems within wide binaries. 
Figure 2a predicts that many should undergo additional dynamical instabilities triggered by 
their stellar companions. Thus, these systems should experience an even greater number of 
planet-planet scattering events than isolated planetary systems. 
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This raises the possibihty that the eccentricities of exoplanets may hold a signature of 
the dynamical process illustrated in Figure 1. Indeed, the overall distribution of exoplanet 
eccentricities provides compelling evidence of our disruptive mechanism. Figure 3a compares 
the observed eccentricity distribution of all Jovian-mass (msinz > 1 Mjup) exoplanets found 
in binarie^ with the distribution of Jovian-mass planets around isolated stars. As can be 
seen, the distribution of planets within wide binaries is significantly hotter than planetary 
systems without known stellar companions. A Komolgorov-Smirnov test returns a probabil- 
ity (or p-value) of only 0.6% that such a poor match between the two datasets will occur if 
they sample the same underlying distribution. Thus, we reject the null hypothesis that the 
distributions are the same. Although it consists of just 20 planets, our wide binary planetary 
sample contains the two most eccentric known exoplanet orbits, HD 80606b and HD 20782b 
(see Figure 2). Furthermore, these excited eccentricities seem to be confined to only very 
wide binary systems. Figure 3a also shows the eccentricity distribution of planets residing 
in binaries with average separations below 10'^ AU. Unlike wider binaries, here we see that 
these eccentricities match very closely with the isolated distribution. (A K-S test returns 
a p- value of 91%.) This suggests the variable nature of distant binary orbits is crucial to 
exciting planetary orbits. Large eccentricities of planets within binaries have previously been 
explained with the Kozai resonanc^^'^'^'^, yet this effect should be most evident in these 
tighter binary systems. 

We perform additional simulations attempting to explain the observed eccentricity exci- 
tation in Figure 3a with the mechanism illustrated in Figure 1. These additional simulation 
sets are summarized in Table 1 (B1-B3). Unlike the internally stable planetary systems in 
simulation set A, these simulated systems consist of 3 approximately Jovian-mass planets 
started in unstable configurations (to induce planet-planet scattering) and evolved for 10 
Gyrs (see supplementary information). In the simulation sets presented in Figure 3b, we 
naturally reproduce both observed eccentricity distributions using the same initial planetary 
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systems. When our planetary systems are run in isolation (set Bl in Table 1) planet-planet 
scattering caused by internal instabilities yields the observed planetary eccentricities for iso- 
lated stars (K-S test p- value of 0.42). Then when a 0.4 binary companion is added to 
each system (set B2 in Table 1) the eccentricity distribution is heated further, and again the 
match to observations is quite good, with a K-S test p-value of 0.46. 

In Figure 3c, the match to observed planetary eccentricities is much poorer. Here we 
rerun our binary simulations with galactic perturbations shut off to yield static binary orbits 
(set B3 in Table 1). In this case, the eccentricity distribution is barely more excited than the 
isolated cases, indicating that the variable nature of wide binary orbits is crucial to heating 
planetary eccentricities. Otherwise, most stellar companions always remain far from the 
planets. 

In Figure 3d we reexamine simulation set B2 to determine which types of planetary 
systems are most influenced by wide binary companions. By examining the planetary systems 
after only 10 Myrs, we can view them after most have experienced internal instabilities but 
before the binary has played a large role (since its effects are delayed). We find that 70% of 
our planetary systems have collapsed to two planets. (The remaining are comprised of nearly 
equal numbers of 1- and 3-planet systems.) We then spht these two-planet systems into those 
with the outer planet beyond 10 AU and those with all planets confined inside 10 AU. In 
Figure 3d, the final {t — 10 Gyrs) eccentricity distribution is shown for both subgroups of 
planetary systems. As can be seen, the more extended planetary systems eventually yield 
much more excited eccentricities compared to the compact systems. This is because binaries 
do not have to evolve to such low pericenters to disrupt extended systems. In fact, the 
observed wide binary planetary eccentricity distribution cannot be matched without using 
wide binaries with planets beyond 10 AU (p = 0.016 from a K-S test). Assuming planets 
form similarly in wide binaries and isolated systems, the planetary eccentricity excitation 
observed within wide binaries may offer new constraints on the bulk properties of isolated 
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giant exoplanet systems, which dominate the giant exoplanet catalog. While most detection 
efforts are currently insensitive to planets with periods beyond ~10 years, our work argues 
that massive longer period planets (beyond ~10 AU) should be common around isolated 
stars. Indeed, such distant planets have recently been directly observed^ and microlensing 
results suggest many such planets reside far from host stard^. 

Due to the variable nature of their orbits, very distant binary companions may affect plan- 
etary evolution at least as strongly as their tighter counterparts. This represents a paradigm 
shift in our understanding of planet-hosting binaries, since previous works tend to assume 
only tighter binaries strongly influence planetary system evolutiorP^I^. Intriguingly, the ec- 
centricities of planets in wide binaries may provide new constraints on the intrinsic architec- 
tures of all planetary systems. To further develop this prospect, searches for common proper 
motion companions to planet-hosting stars should be continued and expandecP^ 
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Name 


Number 


Planet 


Planet 


Binary 


Binary 


External 




of Planets 


Masses 


a-rangc 


Mass 


a* 


Perturbations 






(Mj,p) 


(AU) 


(Mo) 


(AU) 


Included 


A 


4 


SS 


SS 


0.1 - 1.0 


1,000 - 30,000 


Tide + Stars 


Bl 


3 


0.5 - ~15 


2 - ~15 


None 


None 


None 


B2 


3 


0.5 - ~15 


2 - ~15 


0.4 


1,000 - 30,000 


Tide + Stars 


B3 


3 


0.5 - ~15 


2 - -15 


0.4 


1,000 - 30,000 


None 



Table 1 - Initial Conditions of Simulation Sets. SS refers to planetary systems resem- 
bling the solar system's four giant planets, and a represents semimajor axis. "Tide + Stars" 
refers to perturbations from the Galactic tide and passing field stars. 
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Figure 1: Simulation of a binary-triggered instability in our solar system. The pericenter 
and apocenter are plotted for Jupiter (red), Saturn (gold), Uranus (cyan), Neptune (blue). 
The binary's semimajor axis (dotted black line) and pericenter (solid black line) are also 
shown. 
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Figure 2: Disruption as a function of binary mass and separation, a. Map of the fraction of 
systems in set A (see Table 1) that lost at least one planet via instability. Binary mass is plotted on the 
y-axis, while the x-axis marks binary semimajor axis. b. The median binary pericenter below which an 
instability is induced in planetary systems as a function of binary mass and semimajor axis. c. The fraction of 
instabilities that occur after the first 100 Myrs of evolution as a function of binary mass and semimajor axis. 
In each panel black data points mark the masses and presumptive semimajor axes of the HD 80606 and HD 
20782 binaries, which host the two most eccentric known planetary orbitd^Sl 123^ While HD 80606b has been 
reproduced with a Kozai-driven mechanism, this process is markedly slower in even wider binaries such as 
HD 20782lPl. Moreover, the presence of more than one planet suppresses these Kozai oscillations^ ^£91 
However, our disruptive mechanism naturally collapses many systems to one planet, still enabling Kozai 
resonances to contribute to eccentricity excitation. Panel a of this figure suggests that binary-triggered 
instability rates become extremely high as binary semimajor axes drop below ^10'^ AU, which could mean 
that tighter binaries trigger planetary system instabilities even more efficiently than those plotted here. 
However, the initial conditions assumed for both our planetary orbits (solar system-like) and binary orbits 
(isotropic) become questionable for binary semimajorl&es below ^10'^ AU (see supplementary information). 
Another interesting aspect not immediately obvious in panel c is that instability times decreases at the 
largest binary semimajor axes. This is because such binaries are rapidly unbound (or "ionized") by stellar 
impulses, making it impossible for these binaries to trigger instabilities at very late epochs (see supplementary 
information) . 




Figure 3: Eccentricity excitation of planets of wide binaries. Comparison of plane- 
tary eccentricity distributions. Red lines correspond to systems with very wide (a^, > 10^ 
AU) companions and blue are isolated systems. Additionally, solid lines mark observed 
distributions, while dashed lines describe the innermost planets in our simulated systems 
(sets B1-B3 in Table 1). In all distributions, planets with a < 0.1 AU are excluded to 
remove tidally circularized orbits, a. A comparison of observed exoplanet eccentricities 
within tighter (a^, < 10^ AU) binaries (black line) to those observed in very wide (a^, > 10^ 
AU) binaries and isolated systems, b. Eccentricities of simulated 3-planet systems after 
10 Gyrs of evolution (sets Bl and B2 in Table 1). c. Simulations from panel B are rerun 
with no galactic perturbations (set B3 in Table 1). d. The final eccentricities of two dif- 
ferent subgroups of panel b's wide binary simulations: systems that consisted of 2 planets 
extending beyond 10 AU at 10 Myrs (black), and 2-planet systems confined inside 10 AU 
at 10 Myrs (green). One notices in panel b that the presence of a wide binary does not 
seem to enhance the production of very extreme planetary eccentricities. However, 1/4 of 
our wide binary systems have planets driven into the central star (1.7 times the rate within 
isolated systems). Tidal dissipation not included in our models could strand these planets 
in very eccentric orbits before they collide with the central stai'^'^. Interestingly, binaries 
also completely strip 20% of our systems of planets, yielding naked stars that once hosted 
planets (see supplementary information). 
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Supplementary Information 



1 Outline 

Here we present additional details on the work discussed 
in our main paper. This supplementary information is 
organized in the following sections: 

• Section 2 - Simulations Description: Details of 
the numerical algorithms used to perform our sim- 
ulations. This discussion includes our assumed or- 
bital distribution of wide binaries and giant planets 
as well as our assumed initial planetary semimajor 
axes and spacing. 

• Section 3 - Supplementary Discussion: Sev- 
eral abbreviated discussions from the main text 
conceming key results are expounded upon here. 
These include: 

- 3.1: A justification of our choice of initial 
planetary orbits 

- 3.2: How the initial distribution of planetary 
semimajor axes evolves over time. 

- 3.3: Our choice of simulation integration 
times and how our results vary with shorter 
integrations. 

- 3.4: Discussion of biases and contamination 
in the observed exoplanet distributions used 
in the main paper 

- 3.5: A description of the evolution of our iso- 
lated planetary systems and their rate of inter- 
nal instabilities 

- 3.6: A systematic study of how binary- 
triggered instability rates vary with the radial 
distribution of planets 

- 3.7: The fates of planets lost from our simu- 
lations 

- 3.8: Stars that are completely stripped of 
planets by wide binary companions 

- 3.9: The rate that wide binaries are ionized, 
or broken up by stellar pertmbations 



- 3.10: The potential effects of a more realistic, 

time-varying galactic environment 

- 3.11: Exploratory simulations and a discus- 
sion of planetary system evolution within 
tighter binaries 

- 3.12: Discussion of whether the excited plan- 
etary eccentricities in wide binaries could in- 
stead be explained by the wide binary star for- 
mation mechanism 

• Section 4 - Supplementary Table: A list of the 
observed exoplanet systems in wide binaries (and 
their key characteristics) that were used to con- 
struct the observed eccentricity distributions of 
planets within binaries in the main paper. 

2 Simulations Description 

We have run two distinct types of simulations: those 
modeling the effects of a wide binary companion on 
planetary systems like our own solar system (set A in the 
main paper) and those modeling a wide binary's effects 
on generically generated planetary systems (sets B1-B3 
in the main paper). These two batches are described in 
the following two subsections. 

2.1 Solar System Simulations 

For our simulations modehng the evolution of our solar 
system's four giant planets in the presence of a wide bi- 
nary, we used the wide binary integrator included in the 
MERCURY simulation package^'- The giant plan- 
ets were initially placed on their current orbits around 
the Sun in each simulation, while the binary compan- 
ion's mass and initial orbit was varied. Four different bi- 
nary companion masses were modeled: 0.1, 0.3, 0.6 and 
1.0 M0. In addition, thirteen distinct initial semimajor 
axes were used for the binary orbit: 1,000, 1,500,2,000, 
3,000, 4,000, 5,000, 7,500, 10,000, 12,500, 15,000, 
17,500, 20,000, and 30,000 AU. All other initial binary 
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orbital elements were drawn from an isotropic distribu- 
tion. For each mass-semimajor axis combination, 50 
different simulations were run, varying the initial binary 
eccentricity, inclination, argument of pericenter and lon- 
gitude of ascending node. This yielded a total of 2600 
nms. Each simulation was integrated for 10 Gyrs using a 
200-day timestep. Even though the age of the solar sys- 
tem is only 4.6 Gyrs, the Milky Way thin disk is thought 
tobe~10Gyrsold33'34. 

To model galactic perturbations on the orbits of bina- 
ries, we use conditions that mimic the present day solar 
neighborhood. Although binaries near the Sun may have 
explored a large range of galactic environments-'^, we 
choose to hold the galactic environments fixed to limit 
the number of parameters we explore. For the Galac- 
tic tide, we assume a local disk density of 0.1 Mq/pc^, 
and our tidal model contains both a vertical and ra- 
dial component^^. However, the vertical component is 
roughly an order of magnitude stronger than the radial 
one. 

To model the effects of passing stars on binary mem- 
bers, we employ the impulse approximation rather than 
directly integrating their passages^^. This approxima- 
tion is vaUd as long as the orbital period of the binary 
is much longer than the encounter timescale of the pass- 
ing star. This holds for the vast majority of our stel- 
lar encounters. For instance, the shortest binary periods 
in our simulations are ~20,000 years, whereas the en- 
counter timescale of a passing star at 10,000 AU mov- 
ing at 40 km/s is ^1,000 years. (Such an encounter 
would hardly effect our tightest binary orbits, and the 
encounter timescale will be even faster for smaller im- 
pact parameters.) To set the encounter rates of our var- 
ious stellar types (M, O, white dwarfs, etc.) we use 
a local observed mass function^^ combined with local 
velocity data for each stellar category^'. Our stellar 
model assumes the local density of main sequence stars 
is 0.034 Mq/pc^, and a mean dispersion of ^42 km/s. 

While mercury's wide binary integrator is excel- 
lent for modeling perturbations from a distant star as 
well as close encounters between planets, it cannot han- 
dle close encounters between a binary companion and a 
planet. This is not an issue for this set of simulations, 
however, for when the star comes close enough to our 
planetary system for star-planet encounters to arise, the 
planets will be destabilized already. 



2.2 Generic Planetary System Simulations 

In our second set of simulations, we model the evolu- 
tion of artificially generated planetary systems orbited 
by a wide binary companion. To construct these plane- 
tary systems, we first choose a minimxun planetary mass 
for each system that is between 0.5 and 5 Mjup. To 
choose this mass, we randomly draw from a M~^-^ 
distribution"'^. Once a minimum mass is chosen, we 
next assign masses to the other two planets. To do this, 
we use the distribution of the observed mass ratios in 
known multi-planet systems where each planet is more 
massive than 0.5 Mjup. However, we clip this mass ra- 
tio distribution at 8, since larger ratios could potentially 
yield stellar-mass planets. Both the initial and final dis- 
tributions of planetary masses are shown in Figure Sla. 
(Final distributions are used from the systems without a 
wide binary.) In addition, the distribution of mass ratios 
is shown in Figure Sib. As can be seen, the distribu- 
tions do not evolve greatly. There is a small bias to- 
ward higher masses and mass ratios after systems have 
gone unstable. This is to be expected from planet-planet 
scattering"*'. 

Next, we assign orbits to each planet. First, a mini- 
mum semimajor axis is chosen randomly between 2 and 
5 AU. The next planet's semimajor axis is set to a ran- 
dom value between 3.5^.0 mutual Hill radii beyond the 
interior planet. The semimajor axis of the outermost 
planet is then set beyond the middle planet in an analo- 
gous manner. The mass order of the planets is random- 
ized so that no correlation between mass and semima- 
jor axis exists. With planetary mass and semimajor axis 
chosen, we now assign the rest of the orbital elements. 
Eccentricities are randomly drawn from a uniform dis- 
tribution between and 0.01, while inclinations are ran- 
domly drawn from a uniform distribution between and 
0.01 radians. Finally all other orbital elements are ran- 
domly drawn from a uniform distribution between and 
27r. 200 planetary systems were constructed in this man- 
ner. 

Finally, we assign wide binary companions to our 
planetary systems. We set all binary companion masses 
to 0.4 M0, comparable to the typical mass of a star. 
To generate our binary orbital semimajor axes, we ran- 
domly select them from a distribution uniformly dis- 
tributed in log-space between 1,000 and 30,000 AU'*^. 
All other orbital elements are drawn from an isotropic 
distribution. 200 different binary orbits are generated in 
this manner. 

To integrate these systems we use a modified version 
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planet mass ratios 

Figure SI. a. Cumulative distribution of planetary 
masses used in our planetary systems. Blue shows the 
initial distribution, and red shows the final distribution 
surviving after 10 Gyrs of evolution, b. Cumulative 
distribution of planetary mass ratios used in our plan- 
etary systems. Blue shows the initial distribution, and 
red shows the final distribution surviving after 10 Gyrs 
of evolution. The black is the observed distribution of 
mass ratios in systems containing multiple planets with 
msini > 0.5 Mj^p. 

of MERCURY'S wide binary integrator^'- Unlike the 
solar system simulations, we are very interested in the 
final architecture of these systems, so it is essential to 
accurately integrate these systems for the entire 10 Gyrs. 
There are two possible sources of inaccuracies. The first 
occurs when a planet is driven to very low pericenter, at 
which point integrations in democratic heliocentric co- 
ordinates are known to breakdown**"*. To guard against 
this, we set the initial timestep of our simulations to 1/60 
of the innermost planet's orbital period (rather than the 
standard 1/20). Even with this small timestep, extreme 
eccentricities can cause rapid increases in numerical er- 



ror. Therefore, if the innermost planet's eccentricity ex- 
ceeds 0.7, we stop our symplectic integration and switch 
to a Bulirsch-Stoer integrator. This integrator is slower, 
but can accurately integrate very eccentric orbits. As 
mentioned in the previous section another source of 
potential error occurs during close encounters between 
planets and the binary companion. To guard against this, 
we also switch to a Bulirsch-Stoer integrator if the bi- 
nary's pericenter comes within 100 Hill radii of the out- 
ermost planet's semimajor axis (typically of order 100 
AU for our extended systems). Although there may be 
better integrator-switching criteria, this routine consis- 
tently triggers a Bulirsch-Stoer integration well before 
the symplectic integrator breaks down due to star-planet 
approaches or eccentric orbits. 

It is well-known, however, that switching between 
integration schemes many times can also degrade the 
quality of an integration*^. To prevent rapid switches 
between integration methods, once the Bulirsh-Stoer in- 
tegrator is activated we wait 5 Myrs to determine if there 
are favorable conditions to resume a symplectic integra- 
tion (low planetary eccentricities and large binary peri- 
center). If this is the case, we once again resume the 
sympectic integration using a timestep that is 1/60 of 
the innermost planet's current orbital period. Finally, 
we also cap the number of times that our algorithm can 
switch integration routines at 50. Depending on the ini- 
tial conditions (isolated vs. binary star), this limit of 50 
integration switches is only reached in 1-6% of runs. 

To measure the effect of a binary star, each planetary 
simulation is run twice: once in isolation, and once sur- 
rounded by its binary companion. Like the solar system 
simulations, each run is integrated for 10 Gyrs. 

With this modified integration routine, the numerical 
angular momentum error ( ) was held below 5x10^^ 
for all simulations. Energy conservation was poorer. 
■^^^ < 5 X 10^^ has been used as a benchmark in pre- 
vious scattering simulations"*^. It was found that ^10% 
of our simulations did not initially meet this criterion. 
Ultimately, we found close encounters between planets 
to be the main source of energy error. A Bulirsch-Stoer 
integration is employed during these encounters, and the 
error tolerance of these integrations was initially set to 
1 part in 10^'. When we reran our simulations with 
a Bulirsch-Stoer tolerance of 1 part in 10*^ only ^2% 
of our simulations had energy conservation poorer than 
5 X 10"'', with the worst cases being near — 10^^- 
While this is still not ideal, such a small fraction is un- 
likely to alter the global eccentricity distributions of the 
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main paper in any significant manner. 

3 Supplementary Discussion 

3.1 Initial Conditions of Generic Systems 

As mentioned in the main text, the prescription we use 
for generating our generic planetary systems is likely 
overly simplistic. In particular, we have implied that 
all planetary systems form with exactly 3 Jovian-mass 
planets between 2 and -^10 AU in a very unstable con- 
figuration. However, the point of these simulations is 
not to reproduce the exact dynamical evolution of plan- 
etary systems. Instead, they serve as a proof-of-concept 
demonstrating how the presence of a wide binary can 
alter planetary eccentricities. Although the results pre- 
sented in the main paper are likely not unique to just the 
small comer of parameter space our initial conditions 
explore, we nevertheless believe it provides compelling 
support of our mechanism's importance. We have a dy- 
namical model that reproduces the observed exoplanet 
eccentricity distribution around isolated stars. Depend- 
ing on how radially extended our planetary systems are, 
we show the presence of a wide binary companion can 
alter this same dynamical model to yield an acceptable 
match to the observed exoplanets of known wide binary 
systems. 

One other major assumption in oiu" simulations is that 
we implicitly assume that planetary system formation is 
unaffected by the presence of a distant binary, since we 
integrate the same systems with and without a binary 
companion. While it is not certain if this assumption is 
correct, recent observations of protoplanetary disk dust 
emission around binary and isolated stars indicate that 
the disk emission of wide binaries and isolated stars are 
very similar"*^. Furthermore, the median pericenter of 
our initial distant binary orbits is ~1600 AU, well be- 
yond the observed edges of typical circumstellar disks, 
which have typical radii of ^100 AU^^. Indeed it has 
been demonstrated that much closer binaries are un- 
likely to disrupt protoplanetary disks''^. 

3.1.1 Alternative Initial Conditions 

It is important to mention that there are other known 
ways to excite the eccentricities of planets, which are 
affected by our choice of initial conditions. First, if our 
planetary systems are unstable and contain more than 
3 Jovian planets we expect an even greater number of 



planet-planet scattering events to occur. Consequently, 
the surviving planets will have undergone many more 
close encounters with planets that are ultimately ejected, 
and we expect their final orbits to be more eccentric on 
average^". To verify this, we have also integrated 200 
5 -planet systems in isolation using the same mass dis- 
tribution and spacing used in our 3-planet simulations. 
The eccentricity distribution of the innermost surviving 
planet in these simulations is shown in Figure S2a. We 
see that it is markedly more excited than the observed 
eccentricity distribution of planets around isolated stars. 
In fact, it's median eccentricity (0.44) is even higher 
than the median eccentricity of planets within wide bi- 
naries (0.39). Thus, inclusion of additional planets in 
systems within wide binaries could be an alternative ex- 
planation for the observed excitation of these systems. 

In addition, another obvious way to alter the eccen- 
tricity excitation of planetary systems is by varying the 
inherent level of internal dynamical stability. In our 
3-planet simulations presented in the main paper, we 
purposely place planets close to one another to trigger 
planet-planet scattering episodes almost inmiediately. 
This occurs in nearly every system (and our main pa- 
per's simulation analysis ignores the few systems that 
do not go unstable). If our planets were more widely- 
spaced we would expect a smaller fraction of systems 
to undergo planet-planet scattering episodes. Thus, the 
final eccentricity distribution of all systems would be 
markedly cooler. To demonstrate this, we have inte- 
grated integrated an additional set of 200 3-planet sys- 
tems in isolation. These systems have the same mass 
distribution as the 3-planet simulations used in the main 
paper. However, in this new set of simulations the plan- 
ets are separated by 4.0^.5 mutual Hill radii, as op- 
posed to the 3.5^.0 mutual Hill radii separations used 
in the main paper's simulations. The final eccentric- 
ity distribution of these more widely space systems is 
shown in Figure S2b. This is compared with the more 
packed isolated systems used in the main paper (set B 1). 
We see that using more widely space planets drops the 
median eccentricity from 0.26 to 0.16. This is because 
only 2/3 of our more widely spaced systems undergo 
planet-planet scattering episodes, as opposed to 91% of 
the isolated 3-planet systems used in the main paper. 
Therefore, one more way to explain the heightened ec- 
centricities observed among planets of wide binaries is 
that isolated planetary systems form in inherently more 
stable configurations. 

The problem with both of the above explanations for 
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Figure S2. a. Cumulative eccentricity distribution for 
the innermost planets after 10 Gyrs of evolution in sim- 
ulations beginning with five planets on unstable orbits 
(dashed line). Cumulative distribution of observed ex- 
oplanet eccentricities around isolated stars (solid line). 
Only planets with a > 0.1 AU and msini > 1 Mjup 
are included in the observed distribution, b. Final cumu- 
lative eccentricity distributions for the innermost planets 
in two different sets of 3-planet simulations. The plan- 
ets are initially spaced by 4.0^.5 mutual Hill radii in 
the first set (black) and by 3.5-4.0 mutual Hill radii in 
the second set (blue). 

the observed excited eccentricities of planets in wide bi- 
naries is that they imply the process of planet formation 
is significantly different within very wide binaries com- 
pared to isolated stars. However, we show in our work 
that the members of very wide binaries spend most of 
their lives very far from one another. Phases where a bi- 
nary orbit becomes very eccentric tend to be quite short- 
lived. Thus, we find it unlikely that a distant binary com- 
panion has a strong enough impact on planet formation 
to make the initial architecture of these planetary sys- 



tems systematically different from isolated systems. In- 
stead, we think that delayed instabilities caused by the 
time-varying nature of very wide binary orbits is a much 
more plausible explanation. 

3.2 Planetary Semimajor Axes 

We show the initial and final semimajor axis distribu- 
tions in Figure S3. Here we see that the bulk of our 
planets are initially located between 2 and ^10 AU. Af- 
ter internal instabilities occur, this distribution spreads 
out to both smaller and larger semimajor axes. The sur- 
viving inner planets of our simulations tend to perform 
most of the work ejecting and scattering other planets, 
and conservation of energy requires them to move to 
smaller semimajor axes. In addition, other planets can 
be scattered to large semimajor axes without actually 
being ejected, which populates the high-a tail of our dis- 
tribution. 

It should be noted that the large majority of our sim- 
ulated planets still orbit beyond the orbital distance of 
most known Jovian-mass planets (the observed median 
semimajor axis is ^-l.SS AU). We have intentionally 
done this to avoid prohibitive computing time for these 
simulations. However, when constructing the eccentric- 
ity distributions shown in Figure 3 of the main paper, 
we use only the innermost planet of each systems, as the 
outer planets often have decades-long orbital periods. If 
more distant planets reside in known exoplanet systems, 
it is likely that most have not yet been detected by radial 
velocity searches due to their longer orbital periods. In 
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Figure S3. Cumulative distribution of planetary semi- 
major axes used in our planetary systems. Blue shows 
the initial distribution, and red shows the final distribu- 
tion surviving after 10 Gyrs of evolution in isolation. 
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fact, even our innermost simulated planets tend to have 
more distant orbits than observed planets. However, it is 
still valid to compare the two groups' eccentricities, as 
there is currently little evidence that eccentricities gen- 
erated by planet-planet scattering vary much with orbital 
distance**'. 



3.3 Integration Times 

One possible criticism of the work in the main paper is 
our choice to integrate our generic systems for 10 Gyrs. 
As stated above, our set of generic planetary systems is 
designed to be a "proof-of-concept," and the dynami- 
cal footprints of wide binary companions are strongest 
when our planetary systems are oldest. Furthermore, al- 
though determining stellar ages is fraught with uncer- 
tainty, it appears that Milky Way's thin disk began form- 
ing stars ^10 Gyrs ago^^, and that most stars in the solar 
neighborhood are older than the Sun^**. 

Nevertheless, some exoplanet-hosting stars are cer- 
tainly younger than ~10 Gyrs old, so it is important 
to show how our distribution of simulated planetary ec- 
centricities evolves over time. In Figure S4, we show 
how the eccentricity distribution of our planetary sys- 
tems evolves with age. In panel S4a, we plot the me- 
dian eccentricity of the innermost planets of our simu- 
lated systems vs time. As can be seen, little eccentricity 
evolution takes place between 5 and 10 Gyrs. During 
the final 5 Gyrs, the median eccentricity of the inner- 
most planets of our simulated systems only rises from 
0.382 to 0.417. (Note that the median eccentricity of 
observed planets in wide binaries is 0.39.) In compar- 
ison, this value increases from 0.253 to 0.382 between 
t = 10 Myrs and t — 5 Gyrs. In addition, the bottom 
panel of Figure S4 demonstrates that a K-S test com- 
paring our simulated planets with those in wide binaries 
returns only marginally different results during the fi- 
nal 5 Gyrs. Thus, using systems that are only 5 Gyrs 
old does not drastically weaken the agreement between 
our simulated planetary eccentricities and observed ex- 
oplanet systems. Furthermore, one could easily imag- 
ine accelerating the eccentricity evolution seen in Figure 
S4 by using more radially extended planetary systems. 
In addition, eccentricity evolution could be sped up in 
a denser galactic environment where perturbations are 
stronger (see Section 3.10 for a discussion). 
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Figure S4. a. Plot showing the evolution of the median 
eccentricity of the innermost planets in our simulated 3- 
planet systems as a function of time (leftmost data point 
corresponds to t = 10 Myrs). b. Time evolution of 
the p-value of a K-S test when comparing our simulated 
planetary eccentricities with those of observed planets 
in very wide binaries. 

3.4 Further Notes on Observed Eccentric- 
ities 

The catalog of extrasolar planets is far from complete, 
and there exist well-known biases in the observed data 
(detection of short-period, massive planets). It is im- 
portant then to consider if these biases might explain 
the heightened planetary eccentricities seen within wide 
binary star systems, especially since we use this trend 
to support and constrain our numerical work. Because 
many early radial velocity exoplanet searches looked to 
our own solar system as the archetype planetary sys- 
tem, these campaigns originally targeted isolated stars 
and tended to avoid binary systems. 

Most exoplanet-hosting binary stars are located rel- 
atively close to the Sun (this enables the detection of 
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the stellar companions). Consequently, most of these 
exoplanets are discovered with the radial velocity tech- 
nique, since these campaigns tend to target nearby 
stars. Therefore, when comparing planets around iso- 
lated stars to planets within wide binaries, we only con- 
sider planets discovered via radial velocity, as listed on 
exoplanets . org. This also limits our biases to just 
those associated with this detection technique. 

Works cataloging exoplanets within binary star sys- 
tems can be broadly divided into two groups: 1) studies 
that perform follow-up observations on known planet- 
hosting stars to specifically search for unknown stellar 
companions^'' and 2) studies that attempt to 
cross-match the catalog of known exoplanet host stars 
with catalogs of known stellar multiples or common 
proper motion pairs^^' Thus, most exoplanet-hosting 
binaries were originally mistaken for exoplanet-hosting 
isolated stars. Only after exoplanets were discovered 
around these stars were they realized to be members of 
multi-star systems. Although the catalog of exoplanet- 
hosting binaries is subject to observational biases and is 
compiled from a multitude of different works, this cat- 
alog is ultimately a subset of detections from the same 
campaigns that have yielded the more extensive catalog 
of exoplanets around isolated stars with the same detec- 
tion techniques. Therefore, it is not obvious how the 
excited eccentricities of exoplanets within wide binaries 
could arise from an observational bias not present in the 
catalog of exoplanets that orbit isolated stars. 

3.4.1 Eccentricity Discrepancy as a Function of 
Planetary Mass 

In the main text, we only consider observed exoplanets 
with masses above 1 Mj^p. This is done because re- 
cent analyses of exoplanet eccentricities show a well- 
defined division between the eccentricities of planets 
more massive than ^1 Mjup and those less massive than 
1 Mjup^^'"''. While earlier works placed the dividing 
mass at even higher exoplanet values of ^4 Mjup^**- 
our limited catalog of exoplanets within wide binaries 
would quickly be depleted if we increased our mass cut- 
off. However, it is worth considering how the eccentric- 
ity distribution of exoplanets within wide binaries com- 
pares to isolated exoplanet systems when we instead de- 
crease the minimum planetary mass that we consider. 
In the main text, we find a K-S test returns a p-value 
of 0.006 when we use a minimum planetary mass of 1 
Mjup. When we instead consider all planets more mas- 
sive than 0.25 Mj^p the p-value increases to 0.04, and 
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Figure S5. We use a K-S test to compare observed plan- 
etary eccentricities around isolated stars with planetary 
eccentricities found within wide binaries. For this test, 
we only consider planets with masses greater than a 
minimum mass, m„iin- The p-value of a K-S test com- 
paring these two eccentricity distributions is plotted as a 
function of the minimum planetary mass we consider 



the differences between the two exoplanet distributions 
are less severe. 

A plot of how the p-value of a K-S test varies as a 
function of our minimum planetary mass is shown in 
Figure S5. As can be seen, this value increases with 
decreasing planetary mass, climbing above 0.01 when 
we consider planets more massive than 0.6 Mjup and fi- 
nally settling at 0.035 when all planets are considered. 
Although the statistics are poor (we only have 12 ex- 
oplanets whose masses are between and 1 Mjup re- 
siding within wide binaries), this suggests that the dis- 
crepancy between the two eccentricity distributions is 
confined to high-mass planets. This potentially supports 
our explanation that binary-triggered instabilities cause 
this eccentricity excitation. If a planetary system within 
a wide binary happens to avoid a binary-triggered insta- 
bility, then the evolution of planetary orbits should be 
identical to those around isolated stars. On the other 
hand, if the wide binary evolves to a low pericenter and 
greatly excites the planetary orbits, we expect only high- 
mass planets will preferentially survive, while the low- 
mass ones will be preferentially ejected. Thus, the sys- 
tems most harassed by binaries will be largely devoid 
of low-mass planets, while those unaffected by their bi- 
nary companion may host planetary systems similar to 
isolated systems. 
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3.4.2 Catalog Contamination by Undiscovered Bi- 
nary Systems 

As direct imaging studies have demonstrated^^' 
planet-hosting binaries are typically mistaken for iso- 
lated stars before a stellar companion is detected. This 
then raises the prospect that planet-hosting binaries pol- 
lute the catalog of isolated planet-hosting stars and dis- 
tort their planetary orbital distributions. However, there 
is a limit to the level of contamination. As mentioned 
above, many of the original exoplanet searches avoided 
known binaries and targeted stars thought to be iso- 
lated. Consequently, our catalog of isolated planet- 
hosting stars should have a lower frequency of wide bi- 
naries than a random sampling of field stars, and we can 
use the observed rate of wide binaries in the field as 
a maximum estimate of possible contamination. Stud- 
ies using common proper motion pairs and angular 
two-point correlation functions^' have found that bina- 
ries with a* > 10"^ AU account for ^10% of all stars in 
the solar neighborhood. We can then estimate the max- 
imum effect of binary contamination by assuming that 
10% of all isolated planet-hosting are actually binaries. 
To estimate this, we randomly sample the eccentricity 
distribution of exoplanets within wide binaries to gener- 
ate a list of "fake" planetary eccentricities that is 10% as 
extensive as our catalog of planets around isolated stars. 
Adding this list of fake eccentricities to our real catalog 
only changes the median eccentricity of planets around 
isolated stars from 0.25 to ^0.26. Thus, the effect of 
binary contamination seems to be negligible. 

However, there is also the possibility that isolated 
stars once had distant stellar companions that were even- 
tually stripped off by impulses from passing stars. These 
previously binary stars would then be included as iso- 
lated stars. To correct for this effect in our simulated ec- 
centricity distributions of planets within wide binaries, 
we make sure to only analyze simulated planets orbit- 
ing stars that still retain a binary companion. We cannot 
correct for contamination among our observed isolated 
systems, though. Once again, we can estimate the maxi- 
mum effect of this contamination. In our simulations of 
wide binaries, we see that only for a* > AU are the 
majority of binaries ionized over the age of the Milky 
Way's thin disk (see discussion of Figure S8 and Section 
3.9). Further, upon examining Figure 2 of the main pa- 
per, we see that the rate of binary-triggered instabilities 
begins to diminish beyond a* ^ 20, 000 AU because 
these binaries are ionized very rapidly. Thus, the wide 
binaries most likely to affect planets and then get ion- 



ized later have initial semimajor axes between 10,000 
and 20,000 AU. If -10% of stars form with wide bi- 
nary companions, and if wide binary semimajor axes 
follow an Opik distribution, which is distributed uni- 
formly in log-space, then this range of binary semimajor 
axes would correspond to 2-3% of all stars. Therefore, 
by the result of the previous contamination analysis, this 
effect should be negligible as well. 

3.5 Evolution of Isolated Systems 

While the main text focuses on the evolution of our sim- 
ulations that include a binary companion, it is important 
to also document the evolution seen in our isolated 3- 
planet systems. Most of these systems undergo internal 
instabilities very quickly. (Note that we consider an in- 
stability to occur when a planet is lost from the system.) 
In Figure S6, we plot the fraction of systems that have 
gone unstable as a function of time. One sees that the 
large majority (81%) of our systems become unstable 
before t = 10 Myrs. It should also be noted that 9% of 
our isolated planetary systems are stable for the entire 
10-Gyr integration. We assume that planet-planet scat- 
tering is an ubiquitous process in planetary systems, and 
because of this, we do not include these systems in Fig- 
ure S6 or in the eccentricity distributions (both isolated 
and wide binary) in Figure 3 of the main paper. 

In the final part of the main text, we use the plan- 
etary systems that have collapsed to 2-planet systems 
to argue that wide binaries affect planetary systems ex- 
tended beyond 10 AU much more strongly than those 
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Figure S6. Cumulative distribution of the times that iso- 
lated 3-planet systems undergo their first dynamical in- 
stability. Systems that remain stable for 10 Gyrs are ex- 
cluded from this distribution. 
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confined inside 10 AU. To support this claim, we exam- 
ine the instability rate of similar systems in our isolated 
batch of simulations. At t = 10 Myrs, we examine all 
of the 2-planet systems in our isolated runs. For those 
with an outer planet beyond 10 AU, 16% of our sys- 
tems experience another instability at a later time. This 
is a very similar rate to our systems confined inside 10 
AU (19%). For comparison, the instability rate is much 
higher for our extended 2-planet systems that include 
a wide binary. 60% of these systems experience ad- 
ditional instabilities. Meanwhile our compact 2-planet 
systems with wide binaries have instability rates closer 
to the isolated cases, with only 25% of these systems 
experiencing instabilities after t = 10 Myrs. Thus, it 
seems that wide binary companions preferentially desta- 
bilize extended planetary systems. 

Finally, it should be noted that these instability statis- 
tics only apply to the simulations used to construct Fig- 
ure 3 in the main text (simulation sets Bl through B3). 
In Figures 1 and 2 we integrated the four giant planets 
of the solar system in the presence of a very wide bi- 
nary (simulation set A). Obviously, this planetary con- 
figuration is stable indefinitely in the absence of binary 
perturbations. 

3.6 Instability Rate as a Function of Plan- 
etary Semimajor Axis 

In Figure 3 of the main text, we show that planetary sys- 
tems extended beyond 10 AU are much more likely to 
suffer disruption from their binary companion than sys- 
tems with planets all confined within 10 AU. To better 
quantify this trend, we perform an additional set of sim- 
ple simulations with more well-tuned initial conditions. 
These simulations consist of 1,000 runs with 10 differ- 
ent systems of two 1 Mjup planets (100 runs per plane- 
tary system). In each planetary configuration, the plan- 
ets are started on nearly circular, coplanar orbits sepa- 
rated by 6 mutual Hill radii (a stable orbital configura- 
tion). The parameter varied between configurations is 
the semimajor axis of the outer planet, with values rang- 
ing between 2.5 and 100 AU. These systems are then in- 
tegrated in the presence of a 0.4 Mq companion whose 
orbital semimajor axis is drawn from a uniform distri- 
bution in log-space between 1,000 and 30,000 AU. All 
other elements are drawn from isotropic distributions. 

Our systems are then integrated for 10 Gyrs or un- 
til a planet is lost from the system due to a binary- 
triggered instability. The fraction of planetary systems 
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Figure S7. Plot of the fraction of 2-planet systems that 
lose one or more planets due to binary-triggered insta- 
bilities as a function of the initial semimajor axis of the 
outermost planet. 

that have gone unstable is plotted as a function of the 
outer planet's semimajor axis in Figure S7. As can be 
seen, the instability fraction for planetary systems ex- 
tending to 2.5 AU is ^15% but more than doubles by 
10 AU and exceeds 50% by 40 AU. Although this es- 
sentially just measures the rate that the outermost planet 
is excited by a binary, this excitation often cascades in- 
ward in packed planetary systems, exciting the orbits of 
observable planets This process is illustrated in Fig- 
ure 1 of the main paper, where Uranus is actually the 
first planet ejected via binary perturbations even though 
Neptune's orbit is most strongly perturbed. Again, the 
observed eccentricity excitation in known planets of 
wide binaries combined with the rapid rise in instabil- 
ity rates between 2.5 AU and ^10 AU shown in Figure 
S7 suggests that most systems of giant planets extend to 
tens of AU. 

3.7 Lost Planet Fates 

Because we begin our generic planetary systems in un- 
stable configurations, many planets are lost throughout 
the course of our 10-Gyr integrations. In our isolated 
batches, 12% of our lost planets collide with the cen- 
tral star (r < 0.01 AU), 15% suffer collisions with other 
planets, and the remaining 73% are ejected from the sys- 
tem (r > 10^ AU). In our wide binary batch of simu- 
lations, the fractional breakdown is similar, with 14% 
of lost planets colliding with the star, 13% colliding 
with other planets, and 73% ejected. It should also be 
noted that many more planets are lost in our wide binary 
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simulations (60%) vs. our isolated simulations (41%). 
The enhanced planetary loss rate in our simulations with 
wide binaries is why 70% more of these systems expe- 
rience a planetary collision with the central star when 
a wide binary is included, even though the chance of 
collision is still relatively small. This enhanced loss 
rate is not surprising, as previous works have shown 
that close stellar flybys can greatly increase planetary 
ejections^^' In tighter binaries, it has been shown 
that a small fraction of these ejected planets may be ex- 
changed with the other binary member^^ . Unfortunately, 
our simulations may not be capable of measuring this 
statistic, since MERCURY'S wide binary symplectic al- 
gorithm always assumes the planets are in orbit about 
the primary. We intend to study the phenomenon of 
planetary exchange further in the future. 

3.8 Planet Stripping 

As mentioned in the main text, ~20% of our simu- 
lated systems are completely stripped of planets by their 
wide binary companion. In these cases, the binary orbit 
reaches such a low pericenter that even the innermost 
planet is destabilized. In most of our stripping cases, the 
final surviving planet is driven into the star via Kozai os- 
cillations. This occurs for 64% of our stripped systems. 
In the remaining cases the last planet is dynamically un- 
bound from its host star via strong perturbations from 
the close-passing binary companion. Based on these re- 
sults we predict that up to 1 in 5 wide binary stars pos- 
sess no planetary systems at all, even though they may 
have formed planets early in their history. However, as 
discussed in the main text, for the planets driven into 
their stars by Kozai resonances, tidal dissipation during 
pre-coUision close passages between the host star and 
planet may allow some of these planets to survive*"^- 

3.9 Binary Ionization 

As mentioned in the main text, the impact of our most 
distant binaries on planetary evolution is limited by the 
fact that they are quickly ionized by perturbations from 
the local galactic environment. In fact, of our 200 3- 
planet simulations, ^1/4 of our systems lose their bi- 
nary companion after 10 Gyrs of evolution. Because 
these planetary systems would be classified as isolated 
systems after binary ionization, planetary eccentricities 
from these ionized systems are not included in Figure 3 
of the main text. 
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Figure S8. Fraction of binaries that have been ionized 
at t = 1 Gyr (blue), 5 Gyrs (red), and 10 Gyrs (black). 
This fraction is plotted as a function of the initial binary 
semimajor axis. 

Because our solar system simulations integrated large 
numbers of binary stars at specific initial semimajor 
axes, we can use these simulations to study the ioniza- 
tion rate of wide binaries in the local galactic environ- 
ment. In Figure S8, we plot the fraction of binaries that 
become ionized as a function of their initial semimajor 
axis for various times. As can be seen, our population of 
a* = 30, 000 AU binaries is quickly eroded. After 1 Gyr 
of evolution ^ 1/3 of these binaries have been lost, and 
this fraction exceeds 90% after 10 Gyrs. a* ~ 10, 000 
AU is roughly the point at which binary ionization plays 
a large role in the evolution of these systems. Over half 
of systems with initial semimajor axes beyond this value 
are ionized in 10 Gyrs. Inside of this, the ionization rate 
continues to fall, with < 10% of binaries ionized inside 
a* = 4, 000 AU. 

3.10 A Variable Galactic Environment 

As mentioned above, we choose to hold the both the 
population of passing stars and the strength of the Galac- 
tic tide fixed for our simulations. This is done to keep 
our simulations simple and limit the number of parame- 
ters we vary. However, it has been shown recently that 
most stars within spiral galaxies radially migrate over 
large galactocentric distances on Gyr timescales^^- In 
addition, an outward migration appears likely for stars 
with Sun-like kinematics^^. This suggests that many 
wide binaries in the present day solar neighborhood 
once inhabited denser regions of the Galaxy where ex- 
ternal perturbations from passing stars and the Galaxy's 
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tide would have been stronger. This would decrease 
the timescale on which wide binary orbits evolve, and 
consequently, it is likely that wide binary companions 
would have an even greater impact on the planetary sys- 
tem evolution than our simulations suggest in the main 
text. 

As a simple demonstration of this effect, we rerun our 
solar system simulations (set A in the main paper) that 
have a 0. 1 Mq wide binary companion. In the original 
simulations, 35.7% of our planetary systems ejected at 
least one planet. When we double the local galactic disk 
density and stellar density we find that the percentage of 
destabilized planetary systems increases to 44.5%. Fi- 
nally, if we double the local density again this percent- 
age climbs to 47.4%. In a study of the migratory history 
of Sun-like stars, it was found that the mean density of 
the time-integrated local galactic environment is typi- 
cally more than twice as high as the present day solar 
neighborhood^^, supporting the idea that the simulations 
in the main text represent a conservative estimate of the 
importance of wide binary stars on planetary systems. 

3.11 Tighter Binaries 

For our simulations presented in the main paper, we do 
not consider binaries with semimajor axes below lO'^ 
AU. However, in Figure S9 we extend Figure 2 of the 
main paper down to binary semimajor axes of 500 AU. 
Figure S9a indicates that the fraction of planetary sys- 
tems destabilized by binaries becomes extremely high 
(> 90%) in binaries with semimajor axes below lO'^ 
AU. Based on these results, it is even more surprising 
that only planets of very wide (a* > 10"^ AU) are ob- 
served to have excited eccentricities. Such a trend seems 
to contradict this figure. However, Figure S9 may be 
misleading for several reasons. First, in each of the 2800 
simulations that comprise this figure, we use the solar 
system's four giant planets to measure the rate that bi- 
naries trigger planetary instabilities. This assumes that 
the presence of a binary does not affect planet formation 
enough to yield systematically different types of initial 
planetary systems around binaries of different separa- 
tions. Further, by using planetary systems like the so- 
lar system, we are assuming that binaries form initial 
planetary systems like those around isolated stars like 
our Sun. For very wide binaries, this is a reasonable 
assumption, since we show in the main paper that the 
periods in which two binary member stars closely ap- 
proach each other are rare and brief, and therefore un- 
likely to alter the relatively short epoch of planet for- 
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Figure S9. Alternate version of Figure 2 of the main 
paper extended down to binary semimajor axes of 500 
AU. a. Map of the fraction of systems that lost at least 
one planet via instability. Binary mass is plotted on the 
y-axis, while the x-axis marks binary semimajor axis, 
b. The median binary pericenter below which an insta- 
bility is induced in the planetary system as a function 
of binary mass and semimajor axis. c. The fraction of 
all binary-triggered instabilities that occurred after 100 
Myrs as a function of binary mass and binary semima- 
jor axis. In each panel black data points mark the masses 
and presumptive semimajor axes of the HD 80606 and 
HD 20782 binaries, which host the two most eccentric 
known planetary orbits^^' 
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Figure SIO. Fraction of binary-triggered instabilities 
from Figure S9 that are delayed until after t — 100 Myrs 
as a function of binary semimajor axis. 

mation around each star. However, for binaries with 
a* < 10'^ AU this becomes a much more questionable 
assumption. For our binaries with a* — 500 AU, the 
median initial pericenter of their orbits is only ^150 
AU. Thus, they potentially have a much more signifi- 
cant impact on planet formation, and it's not clear that 
we should even expect solar system-like architectures to 
arise in binaries with a* < lO'^ AU. Perhaps then it is 
not very surprising that such planetary systems are al- 
most always immediately destabilized when placed in 
binaries with a* < 10'^ AU. 

The constant significant gravitational interaction be- 
tween members of tight binaries is evidenced in Figure 
S9c. Here we see that the timescale for binaries to trig- 
ger planetary system disruptions becomes much shorter 
below a, of 1,000-2,000 AU. In Figure SIO, we sum 
these simulations over all binary masses to see what 
fraction of instabilities occur well after planet formation 
(t > 10^ yrs) solely as a function of binary semima- 
jor axis. Here we see that for binary semimajor axes 
beyond --3,000 AU at least --95% of all instabilities 
are delayed. Inside a* < lO'^ AU, however, this frac- 
tion plummets to just under 10%, implying that most 
systems are immediately unstable. If these planetary 
systems are instantly unstable, we may actually expect 
planet formation to be altered enough to prevent forma- 
tion of solar system-like planetary architectures in such 
binaries. 

Another reason the low-a* region of Figure S9 may 
be inaccurate is because of the initial distribution of bi- 
nary orbital elements we assume in our simulations. It 
is thought that binaries with a* < lO'^ AU have a funda- 



mentally different formation process than binaries with 
a* > 10^ AU. For binaries with a* > 10^ AU, it is 
largely thought the formation mechanism is dynamical 
capture either due to 3-body interactions or star clus- 
ter dispersal™' Such processes yield an isotropic 
distribution of binary eccentricities and other orbital 
elements^', which is what we assume in our simula- 
tions' initial conditions. On the other hand, tighter bi- 
naries (a* < 10'' AU) are thought to form via fragmen- 
tation of collapsing cloud cores^^. In this scenario, it is 
less clear what the resulting distribution of orbital ele- 
ments should be. Indeed, studies of imaging polarimetry 
of young binaries indicate that the circumstellar disks of 
binary members are coplanar to within 20° for binary 
separations between 200-1,000 AU^'*'^^. This would 
imply that the binary orbital planes must also be nearly 
coplanar with the circumstellar disks, or else these disks 
would rapidly become misaligned. 

Thus, if inclinations between tighter binary orbits and 
planetary disks are colder, the eccentricities may be as 
well. This could also explain why the observed plane- 
tary eccentricities with tighter binaries are colder than 
those found in binaries with a* > 10'^ AU. Such an or- 
bital distribution is much less likely to result in planetary 
system disruptions. To demonstrate this, we integrate 
nearly coplanar (i < 20°), nearly circular (e < 0.1) bi- 
nary companions with masses of 0.4 Mq about the four 
giant planets of the solar system for 10 Gyrs. 

We perform 20 integrations at binary semimajor axes 
of 300, 500, 750, 1,000, 3,000, 5,000, and 10,000 AU 
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Figure Sll. Fraction of systems that lose at least one 
planet due to a binary-triggered instability as a function 
of binary semimajor axis. Systems are run in the pres- 
ence of galactic perturbations {blue) and with perturba- 
tions shut off (green). 
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for a total of 140 runs. The fraction of planetary sys- 
tems that have undergone an instability resulting in at 
least one planetary ejection is plotted as a function of 
binary semimajor axis in Figure Sll. In this plot we 
see that binaries with a* ~ 300 AU are still quite dis- 
ruptive to systems Uke our solar system. However, be- 
tween a* ^ 500 AU and a* ~ 3, 000 AU the plane- 
tary systems are relatively unaffected by the presence of 
the binary. Beyond a* ~ 3,000 AU, the disruption rate 
picks up again, as perturbations from the Galactic tide 
and passing stars are powerful enough to completely 
transform the initial eccentricity distribution. (It should 
be noted that our initial cold eccentricity distribution 
largely suppresses binary-triggered disruptions for bi- 
naries with 1, 000 < a* < 5, 000 AU compared to the 
results of the main paper. This is because the timescale 
for the Galactic tide drive a circular orbit to an extreme 
eccentricity is longer than our integration time. Plane- 
tary system disruption in these binaries typically results 
instead from initially moderately eccentric binary orbits 
driven to extreme eccentricities.) 

Although these simulations' initial conditions are ide- 
alized, they illustrate an important point. There are 
zones of binary orbital parameter space that allow cer- 
tain classes of planetary systems to remain stable indef- 
initely in tighter binaries. When the binary separation is 
increased, however, galactic perturbations drive binaries 
far from their initial orbital configurations, regularly dis- 
rupting the same planetary systems. To further illustrate 
this, we rerun our simulations with galactic perturba- 
tions turned off. The planetary instability rates found in 
these simulations are also shown in Figure Sll. As can 
be seen, without galactic perturbations, the instabilities 
occurring within the largest binary semimajor axes have 
disappeared. 

We admit that our explanations offered for the unex- 
cited planetary eccentricities observed in tighter bina- 
ries are not fully tested and are a bit speculative. How- 
ever, there are already many works focused on plane- 
tary dynamics for binary separations at or below ^10^ 
AU^*' while our work is among the first to focus 

on planets within the widest binaries. Planetary dynam- 
ics within tighter binaries will be a topic for our future 
research, but a fully detailed study of this problem is 
beyond the scope of the present work. 

3.12 Wide Binary Formation 

In the main paper, we demonstrate that the orbital eccen- 
tricities of observed giant exoplanets with wide binary 



companions are statistically higher than those without 
stellar companions. We argue that this is due to the dy- 
namical evolution of the wide binary systems under the 
influence of other passing stars and the Galactic tide. 
However, the exact process of wide binary formation re- 
mains poorly understood, and here we consider the pos- 
sibility that the process of wide binary formation yields 
more excited planetary systems. 

Several recent works have proposed that wide bina- 
ries form during the dissolution of star clusters^'- ™''^. 
In such scenarios, the two stars are already tenuously 
bound during cluster dissolution and are "frozen out," 
or they become bound to one another as they escape the 
dissolving cluster potential. These formation mecha- 
nisms seem unlikely to excite planetary eccentricities in 
these systems because the two member stars should only 
rarely closely encounter one another during the forma- 
tion process. Such encounters are not required in these 
formation scenarios, and they would only arise if a sys- 
tem forms with an initial binary orbit that is extremely 
eccentric. Such initial orbits are unlikely for any plausi- 
ble binary orbit eccentricity distribution, even isotropic. 

However, another new formation mechanism has 
been proposed suggesting that most wide binaries are 
bom as triple star systems^^. In this mechanism, the 
triple system quickly goes unstable. Typically this 
would result in the ejection of one of the stars. How- 
ever, if the instability occurs quickly enough the poten- 
tial of the star-forming molecular cloud core can instead 
trap a star in a distant orbit before it is ejected. What 
results is a "binary" that is actually composed of one 
well-separated star and two other stars in a much lighter 
orbit about each other. Unlike the previous binary for- 
mation mechanisms, close encounters between member 
stars could be conmion as triple star systems go unsta- 
ble in this scenario. Thus, it seems possible that any 
planetary systems that are present could be significantly 
perturbed during these initial triple star system instabil- 
ities. 

Perhaps then, the excited planetary eccentricities we 
discuss in the main paper are actually a signature of 
wide binary formation from triple stars rather than a sig- 
nature of the subsequent dynamical evolution of wide 
binaries. There are two reasons we do not believe this is 
the case, however. First, if most wide binaries form from 
triple systems this requires most known wide binaries to 
actually be comprised of three stars (an isolated star and 
a tight binary). In the main paper, we analyze the orbits 
of 20 known exoplanets with wide binary companions. 
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Of these 20, only 4 reside in systems where the distant 
binary companion is known to be a much tighter binary 
itself. While these are well-studied systems, it is pos- 
sible that the tight binary has been mistaken for a sin- 
gle star in some of the more distant systems. However, 
systems within ~30 pc should indeed have all stellar 
members detected^^. Of the known planet-hosting wide 
binaries, six are actually within ~30 pc, and only one 
of these systems (16 Cygni) is known to be actually be 
comprised of three stars, which suggests that most of 
the more distant planet-hosting wide binaries also only 
have two stellar components. Thus, the predominance 
of conventional "two-star" wide binaries in our plane- 
tary sample suggests that triple stars do not yield most 
wide binaries that host giant planets. (It should also be 
noted that these 4 planets are not particularly eccentric 
members of the overall sample. Only one has an eccen- 
tricity that exceeds the overall sample's median.) 

A second reason we believe that instabilities of triple 
star systems do not explain the heightened eccentricities 
of exoplanets within wide binaries relates to the tim- 
ing of triple star instabilities. Because the potential of 
a molecular cloud core is necessary to prevent the dis- 
tant star from being ejected, these instabilities must oc- 
cur very early in the lifetime of stars, and indeed most 
are predicted to occur in the protostellar stage^^. Con- 
sequently, it is not at all clear that the stars involved in 
this process have formed planets yet when the instability 
occurs. If planetary systems do not form until after this 
instability then we expect the distant stellar companions 
to have little effect on the final orbital eccentricities of 
the planetary systems. Only later when galactic pertur- 
bations drive the stellar orbit to very low pericenter as 
described in the main paper can the distant companion 
excite or disrupt the planetary system. 



suming an isotropic distribution of orbital elements™. 
This is just 1.26 times larger than the projected sepa- 
ration, and we use this value to divide our sample into 
planets of very wide binaries (a* > 10"^ AU) and plan- 
ets of tight binaries (a* < 10'^ AU). Based on this di- 
vision, we have 20 planets within very wide binaries 
and 23 planets that reside in tight binaries. As stated in 
the main text, our analysis does not include planets with 
msini < 1 Mj^p or a < 0.1 AU. These cuts exclude 
~43% of all known planets within binaries. 



4 Supplementary Table 

In Figure 3 of the main paper, we construct two dis- 
tribution functions of observed exoplanet eccentricities 
within binary star systems. Those planets are listed on 
the next page in Table SI along with their orbits and 
masses as well as a description of their binary star sys- 
tem. (It should be noted that in a few cases, the distant 
binary companion is itself actually a much tighter binary 
system. For large separations, the potential of this bi- 
nary will closely resemble a point mass.) Although only 
the projected separations of most binaries are known, 
the most probable semimajor axis can be calculated as- 
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Name 


Projected 


a* 


msini 


a 


e 




Stellar Separation 












(AU) 


(AU) 


(Mjup) 


(AU) 




HD 38529c 


12,040 


15,180 


17.7 


3.695 


0.36 


HD 20782b 


9,080 


11,550 


1.8 


1.36 


0.92 


HD 40979 Ab 


6,400 


8,070 


3.83 


0.855 


0.269 


HD 147513b 


5,360 


6,760 


1.21 


1.32 


0.26 


HD 222582b 


4,750 


5,990 


7.75 


1.35 


0.76 


HD 125612b 


4,750 


5,990 


3.0 


1.37 


0.46 


HD 12561 2d 


4,750 


5,990 


7.2 


4.2 


0.28 


HD 213240b 


3,900 


4,920 


4.72 


1.92 


0.421 


GJ 777b 


3,000 


3,780 


1.56 


4.01 


0.313 


HD 7449b 


2,960 


3,730 


1.1 


2.3 


0.82 


HD 7449c 


2,960 


3,730 


2.0 


4.96 


0.53 


HD 89744b 


2,460 


3,100 


8.58 


0.934 


0.677 


HD 219449b 


2,250 


2,840 


2.9 


0.3 


0.0 


HD 80606b 


2,000 


2,520 


3.9 


0.453 


0.9336 


30 Arietis Bb 


1,500 


1,890 


9.88 


0.995 


0.289 


55 Cancri d 


1,065 


1,340 


3.82 


5.74 


0.014 


1 1 Comae b 


995 


1,250 


19.4 


1.29 


0.231 


16 Cygni b 


860 


1,080 


1.68 


1.681 


0.681 


HD 142022b 


820 


1,030 


4.47 


2.93 


0.53 


GJ 676b 


800 


1,010 


4.9 


1.82 


0.326 


HD 178911b 


785 


990 


7.35 


0.345 


0.139 


Upsilon And c 


750 


945 


13.98 


0.832 


0.224 


Upsilon And d 


750 


945 


10.25 


2.53 


0.267 


Upsilon And e 


750 


945 


1.059 


5.24 


0.005 


HD 188015b 


684 


862 


1.5 


1.203 


0.137 


HD 196050b 


510 


643 


2.9 


2.54 


0.228 


HD 132563 Bb 


397 


500 


1.49 


2.6 


0.22 


HD 65216b 


253 


319 


1.21 


1.37 


0.41 


HD 156846b 


250 


315 


10.45 


0.993 


0.8472 


HD 27442b 


240 


302 


1.56 


1.271 


0.06 


HD 28254b 


238 


300 


1.16 


2.15 


0.81 


GJ 667 Cc 


204 


257 


4.54 


0.123 


0.27 


Gamma Leo b 


170 


214 


8.78 


1.19 


0.14 


HD 19994b 


151 


190 


1.68 


1.42 


0.3 


HD 195019b 


150 


189 


3.69 


0.1388 


0.0138 


HD 142b 


138 


174 


1.31 


1.045 


0.26 


HD 1 14762b 


130 


164 


11.68 


0.363 


0.3359 


HD 177830b 


97 


122 


1.49 


1.22 


0.009 


GJ 3021b 


68 


86 


3.37 


0.49 


0.511 


HD 196885b 


24 


30 


2.96 


2.6 


0.48 


G186b 


21 


26 


3.91 


0.113 


0.0416 


HD 41004 Ab 


20 


25 


2.56 


1.7 


0.74 


Gamma Cep Ab 


20 


25 


1.6 


2.044 


0.115 



Table S 1 . List of the planets residing in binaries which are used to construct the observed distribution functions in 
Figure 3 of the main paper Columns from left to right are: planet/star name, binary projected separation distance, 
most likely binary semimajor axis, planetary mass, planetary semimajor axis, and planetary eccentricity. Systems 
are listed in order of decreasing binary separation. 
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